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The effect of a pressure gradient on the local heating disturbance of rectangular cavities 
tested at hypersonic freestream conditions has been globally assessed using the two-color 
phosphor thermography method. These experiments were conducted in the Langley 31-Inch 
Mach 10 Tunnel and were initiated in support of the Space Shuttle Return-To-Flight 
Program. Two blunted-nose test surface geometries were developed, including an expansion 
plate test surface with nearly constant negative pressure gradient and a flat plate surface 
with nearly zero pressure gradient. The test surface designs and flow characterizations were 
performed using two-dimensional laminar computational methods, while the experimental 
boundary layer state conditions were inferred using the measured heating distributions. 
Three-dimensional computational predictions of the entire model geometry were used as a 
check on the design process. Both open-flow and closed-flow cavities were tested on each test 
surface. The cavity design parameters and the test condition matrix were established using 
the computational predictions. Preliminary conclusions based on an analysis of only the 
cavity centerline data indicate that the presence of the pressure gradient did not alter the 
open cavity heating for laminar-entry/laminar-exit flows, but did raise the average floor 
heating for closed cavities. The results of these risk-reduction studies will be used to 
formulate a heating assessment of potential damage scenarios occurring during future Space 
Shuttle flights. 


I. Nomenclature 


BE 

Bump Factor 

BE 

LJ1 max 

maximum Bump Factor on the cavity end wall 

c P 

pressure coefficient 

h 

heat transfer coefficient, h=q/(H aw -H w ), (lbm/ft 1 2 /s) 

H 

enthalpy (btu/lbm) 

L, W, H 

cavity length, width, and depth (in) 

M 

Mach number 

P 

pressure (psi) 

q 

surface heat transfer rate (btu/ft 2 /s) 

Rn 

model reference nose radius (in) 

Re 

unit Reynolds number (1/ft) 

Re e 

momentum thickness Reynolds number 
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RTF Return-to -Flight 

s distance along the streamline (ft) 

T temperature (°R) 

t time (s) 

TGP thermographic phosphor 

U velocity magnitude (ft/s) 

x axial distance from model leading edge (in) 

y spanwise distance from model centerline (in) 

z distance normal to x-y plane (in) 

a angle- of- attack (deg) 

y ratio of specific heats 

8 boundary layer thickness (in) 

6 boundary layer momentum thickness (in) 

Q Normalized temperature in heating equation 

A Normalized time in heating equation 

p density (lbm/ft 2 3 ) 

Subscripts 

oo freestream static conditions 

avg average 

aw adiabatic wall conditions 

FR conditions from Fay-Riddell calculation for a hemisphere 

local local condition 

ref reference condition 

T1 reservoir conditions 

T2 stagnation conditions behind a normal shock 

w wall conditions 


II. Introduction 

The Final Report of the Columbia Accident Investigation Board (CAIB) was released in August 2003. The CAIB 
identified “a breach in the Thermal Protection System of the leading edge of the left wing, caused by a piece of 
insulating foam ...” as the probable event resulting in the loss of the Space Shuttle Columbia during flight STS-107 
on February 1, 2003. Many possibilities were investigated prior to reaching this assessment, among them the 
creation of imp act- induced cavities in the thermal protection system (TPS) tiles, resulting in local augmentation of 
the heating and eventual burn through of the aluminum substructure. Though impact-induced damage to the tiles 
was eventually excluded as the cause of the Columbia accident, the CAIB recommended that NASA “Develop, 
validate, and maintain physics-based computer models to evaluate Thermal Protection System damage from debris 
impacts. These tools should provide realistic and timely estimates of any impact damage from possible debris from 
any source that may ultimately impact the Orbiter. Establish impact damage thresholds that trigger responsive 
corrective action, such as on-orbit inspection and repair, when indicated .” Langley Research Center, operating in 
concert with the Shuttle Accident Investigation Team and the Return-To-Flight Program, has provided an 
experimental cavity heating database of nearly 775 wind tunnel runs, along with many computational simulations, of 
this complex fluid dynamic and aeroheating environment for analysis and modeling by the Damage Assessment 
Aeroheating Team at Johnson Space Center. The predominance of these data and simulations are for non-surface- 
breaching, impact damage situations where near- zero pressure gradients occur, such as on the windward surface tile 
acreage. Because of the fast pace of this effort, most of these results are only in the earliest stages of reporting and, 
accordingly, are not readily available to the general community at this time. The purpose of this paper is to present 
one aspect of this effort; that is the impact on local heating augmentation of an expanding flow, resulting in a 
negative pressure gradient, over a cavity damage site. Based on the outcome of the analysis, allowances for these 
effects may be required in the engineering modeling tools currently being developed to assess potential heating 
augmentation resulting from damaged surfaces. Results from these analyses directly impact the go/no-go decision 
required prior to any spacewalk tile repair mission performed by the shuttle crew. 
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III. Literature Assessment 

The supersonic/hypersonic cavity flow literature (particularly that for laminar test conditions), while helpful, is 
particularly sparse on both the local and downstream effects of many fundamental flow parameters. Fletcher, et al. 
published the survey paper “A Review of Heat Transfer in Separated and Reattached Flows” in 1970. Nestler 
updated this survey with more recent work in 1985 in “The Effects of Surface Discontinuities on Convective Heat 
Transfer in Hypersonic Flow.” Together, these papers include much of the existing, pertinent cavity flow literature. 
They cited theoretical models of cavity flow developed by Burggraf (1965), Chapman (1956), Carlson (1959), 
Chang (1966), and Lamb (1967); however, none of the theories adequately capture the physics of laminar and 
turbulent three-dimensional cavity flows as will be demonstrated via global surface heating distributions presented 
herein, though they offer a framework for experimental studies. Numerical studies were conducted by Adams 
(1974), Morgenstern and Chokani (1994), Zhang, et al. (2001), and more recently for the CAIB and the Shuttle RTF 
program by Wood, et al. (2004) and Pulsonetti, et al. (2005). The computational work presented by Wood uses the 
previously mentioned unpublished RTF experimental open cavity results for CFD comparison and validation of the 
numerical modeling, while Pulsonetti’ s work deals with flight traceability. Experimentally, laminar two-dimensional 
flows were addressed by Galenter (1975) and Nestler (1966); laminar axisymmetric flows were addressed by Nestler 
(1970), Nicoll (1963, 1964), and Wyborny, et al. (1967); and, laminar three-dimensional flows were examined by 
Cheatwood, et al. (2001), Hahn (1969), and Nestler (1981). Experimental studies with cavities tested in transitional 
boundary layers or to determine the onset of transition were presented by Charbonnier and Boerrigter (1993), 
Boerrigter and Charbonnier (1996), Hollis and Liechty (2002, 2006), Larson and Keating (1960), Liechty, et al. 
(2006a), and Rhudy and Magnan (1964). Turbulent two-dimensional experiments were presented by Charwat, et al. 
(1961a, 1961b), Chin and Seban (1972), Emery (1969), Hunt (1977), Lamb (1968, 1980a, 1980b), Nestler (1968, 
1969, 1970a, 1970b), Shchukin, et al. (1980), and Stallings and Wilcox (1987); turbulent axisymmetric experiments 
were conducted by Hunt (1974) and Netterfield (1989); and turbulent three-dimensional experiments were published 
by Wilcox (1990, 1991). In general, the heating measurements were obtained with sparsely-spaced discrete sensors, 
the exception being those global phosphor thermography measurements presented by Cheatwood, et al., and Hollis 
and Liechty for circular cavities. Because of this spatial measurement sparseness, much of the three-dimensional 
nature of the surface heating profile is missed, as is in many cases the important peak heating value. While most of 
the idealized cavity geometries are rectangular with a flat bottom profile (a few have an arc bottom), the paper by 
Emery is of particular significance because heating profiles are presented for models with notched, cutback profile 
geometries that may be more representative of potential impact damage conditions. Two papers by Blair and 
Stallings (1986) and Stallings, et al. (1991) are significant because they present oil flow visualizations in and around 
the cavity and vapor screen visualizations of the cavity crossflow plane, clearly showing the growth and 
development of the cavity vortex structure. The important class of cavities known as gaps (i.e., cavities with length- 
to-depth ratio less than 1) has been examined by Coats, et al. (1970), Johnson (1973), Throckmorton (1974), and 
Weinstein, et al. (1975). 

As demonstrated above, a large amount of cavity research studies are available in the published literature. 
However, there exists only a limited amount of laminar flow experimental data for evaluating heating effects in 
supersonic/hypersonic cavities, and most of the existing parametric variations and correlations have been developed 
using turbulent heating measurements. Further, turbulent methods have typically been used for thermal assessments 
of damage on the Space Shuttle, even though much of the re-entry trajectory provides laminar edge conditions on 
the vehicle. The use of these turbulent methods may be overly conservative for the laminar conditions. 

No evidence of the effect of pressure gradient on the heating was found in the literature. An objective of the 
present study is to provide sufficient cavity heating data to assess the effects of pressure gradient under laminar 
entry flow conditions. Specifically, the focused objective is to address the binary yes-no question of whether 
pressure gradient is an important flow parameter for cavity damage. A secondary objective is to provide enhanced 
surface heating information under controlled conditions to the CFD community for flow simulation studies. 

IV. Supersonic/Hypersonic Cavity Flow Physics 

Based on the literature survey, the following overview of cavity flow physics is presented. Length-to- depth ratio 
is typically used to distinguish between and classify different cavity flow regimes, as depicted in Figure 1. Very 
short or deep cavities with L/H < 1 are known as gaps. In this case, shearing induced by the main stream flow causes 
the development of a column of counter-rotating vortices within the gap numbering approximately H/L. Alternating 
hot spots are developed in the gap when the vortices directionally align and impinge on the sidewall. Two stable 
flow conditions exist for cavities with length L/H > 1. The first is a short cavity in the length range 1 < L/H <10, 
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known as an open cavity in the literature. The mainstream flow does not enter the cavity directly, though there may 
be some mass interchange with the low-energy vortical flow inside the cavity. In this case, the pressure in the cavity 
is typically above the ambient and climbs to a peak at the downstream lip. The heating drops to values significantly 
below the undisturbed value and rises slowly to a peak value on the downstream lip. The other stable solution is a 
long cavity, also known as a closed cavity, with length L/H > 14. In this case, three distinct flows may develop if the 
cavity is long enough. First, the upstream flow is able to turn into the cavity and impinge on the floor, creating an 
aft-facing-step flow field. Next, a boundary layer on the floor may develop and recover to the ambient level outside 
the cavity. Finally, as the flow approaches the end wall it will turn outward and create a forward-facing-step flow 
field. For long, deep cavities, the pressure gradients may be severe where the flow turns and strong expansion and 
shock waves will be generated. Viscous shearing generated by this flow turning will augment the heating to levels 
significantly higher than the ambient levels on both the cavity floor and the end wall. The pressure in these long 
cavities will decrease below the ambient and steadily increase downstream, reaching large values of over-pressure 
behind the shocks. Vortices will develop on the cavity sidewalls as the flow expands around the corner into the 
cavity and on the floor after flow impingement, further augmenting the heating, which may extend laterally around 
the cavity in the most severe cases. These vortices will interact with the cavity end wall and spill into the 
downstream region beyond the cavity. Given laminar inflow, analysis of the heating profiles indicates that most any 
type of outflow may occur, depending on the cavity and its environmental state. The in-cavity flow may remain 
laminar, become transitional, or transition to fully turbulent flow; the downstream possibilities are equally varied 
and are currently the subject of extensive analysis and testing (Liechty 2006a) in support of Shuttle RTF. Cavities in 
the range 10 < L/H < 14 are typically unsteady as the flow alternates between the two bounding conditions; these are 
known as transitional cavities. Transitional cavities are avoided where possible in the present tests because of the 
complexity of the required instrumentation and test time necessary to address flow steadiness. 

The boundaries between the different cavity flow regimes are nominal, at best. For example, different 
researchers have measured L/H values ranging from 9 to 1 1 as the upper limit for open cavities and from 12 to 15 as 
the lower limit for closed cavity flow. These limits should therefore be taken only as a guide. Also, it is important to 
note that most of the reported cavity data were acquired in air (y^l.4) with some in helium (y =1.67). Since turning 
angle is a function of both Mach number and y, it is conceivable that the open/closed cavity boundaries will vary 
during re-entry as vehicle boundary layer edge conditions change. At present, experimental data are insufficient to 
estimate the strength of this effect on the aeroheating. 

Gap Open Cavity Transitional Cavity Closed Cavity 

L/H < 1 1 < L/H <10 10 < L/H <-14 L/H > -14 


Boundary Shear 



Figure 1.- Cavity flow regimes. 


V. Establishing the Testing Environment 

When few (if any) direct measurements of the cavity geometry are available and when the exact location of the 
cavity (implying local flow conditions) is unknown, defining the heating environment of a damage-scenario cavity 
presents a significant challenge. Many possible variables and parameters exist. Following the literature, the gross 
geometric parameter space includes: length (L), width (W), and depth (H), and planform, cross-sectional, and profile 
variations. Sidewall and upstream end-wall entry angles may play an important role on the flow expansion into the 
cavity, while sidewall and end-wall exit angles may affect recompression on the downstream wall and further 
expansion around the end-wall corner - all parameters having a direct influence on the heating. The depth 
distribution along the cavity profile may potentially affect how or if the flow enters the cavity and impinges on the 
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floor. The impact of roughness caused by irregular surfaces, protuberances, gap fillers between tiles (present or 
missing), embedded objects, etc. in and upstream of the cavity is unknown, but the roughness may significantly 
impact the shear layer/boundary layer transition process (Liechty, et al. 2006b). Local flow conditions and 
parameters include boundary layer thickness (8) and momentum thickness (0), edge Mach number (M e ), momentum 
thickness Reynolds number (Re 0 ), stream direction with respect to the cavity major axis, chemistry effects reflected 
in the ratio of specific heats (y), wall temperature ratio (T w /T aw ), boundary layer state (laminar, transitional, or 
turbulent) entering/leaving the cavity, and, of course, pressure gradient. 

Representative surface flow conditions for the Shuttle Orbiter (Reuther, et al. 2004) have been computed using 
the DPLR code (Wright, et al. 1998) for the STS- 107 re-entry trajectory. Sample conditions from the trajectory 
computations are presented in Table 1 for nominal altitudes of 243,000 ft (Mach 25) and 201,767 ft (Mach 18). 
Three windward surface streamlines spanning the chine region of the wing are presented in Fig. 2a and Fig. 2b, 
pressure coefficients are presented in Fig. 2c and Fig. 2d, pressure gradient distributions are presented in Fig. 2e and 
Fig. 2f, and local Mach number and local momentum thickness Reynolds number are presented in Fig. 2g and Fig. 
2h, and Fig. 2i and Fig. 2j, respectively. The location of the wing leading edge is most readily apparent in the 
pressure gradient as the sharp negative peak that occurs at approximately 36, 49, and 67 feet for streamlines 1, 2, 
and 3, respectively. Immediately aft of the leading edge, a representative dC p /ds value of about -0.15 per foot 
(-0.013 per inch) occurs on the chine. Edge Mach numbers are in the range 2.5 to 2.9. Momentum thickness 
Reynolds number ranges are 80 to 130 for CFD Case 1 and 220 to 340 for CFD Case 6. Therefore, these target 
conditions are selected as representative for designing the experiment. Note, again, that the objective is not to match 
surface conditions exactly, but to assess the effects of pressure gradient at conditions representative of those on the 
surface. 


Table 1.- Selected STS-107 trajectory conditions for CFD modeling. 


CFD 

Case. 

Altitude 

(ft) 

Velocity 

(ft/s) 

Density 

(lb/ft 3 ) 

Temp 

(°R) 

Mach 

Number 

Reynolds 
Number (1/ft) 

AOA 

(deg) 

1 

243,000 

24,116 

2.4350e-6 

391 

24.87 

63189 

40.17 

6 

200,767 

18,164 

1.6293e-5 

429 

17.88 

294950 

39.02 




i i 

a) Case 1 - streamlines. b) Case 6 - streamlines. 

Figure 2.- STS-107 trajectory CFD along streamlines. 


r 
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dC/ds, 1/ft P - psi 




c) Case 1 - pressure. 


d) Case 6 - pressure. 




e) Case 1 - pressure gradient. 


f) Case 6 - pressure gradient. 




h) Case 1 - M e . 


i) Case 6 - M e . 


Figure 2.-Continued. 
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j) Case 1 - Re e . k) Case 6 - Re e , 


Figure 2.-Concluded. 


Given the infinity of potential geometric variations for impact damage, the Shuttle aeroheating community has 
defined a simplified cavity representation. This rectangular geometry (called the “shoebox” geometry due to its 
shape) consists of a flat bottom bounded by straight sides and end walls. Entry (upstream) and exit (downstream) 
and side wall angles are prescribed for the general case. However, for the sake of simplicity, vertical sides and ends 
are used for the present designs. According to existing literature, the governing parameters for a cavity are the 
length- to-depth ratio (L/H), and the depth- to -boundary- layer- thickness ratio (H/8). Lack of detailed width-effect 
information in the literature leads to the decision to scale the width by the depth (W/H). Two different cavity lengths 
were selected, both legacy lengths from previous testing. The first is a short cavity of length L/H=7.2 that is 
representative of an open cavity; and, the second is a long cavity of length L/H=20, representative of a closed cavity. 
Prior testing experience and published information in the literature indicate that the heating is more pronounced in 
closed, deep cavities that allow a significant expansion of the flow into the cavity, particularly when there is 
impingement on the floor. Accordingly, to match existing, unpublished data, a target cavity depth of H/8=l.l is 
selected. Additionally, a larger H/8 translates (via L/H) into a longer cavity, L, enabling better optical access for the 
measurement system. Also, even though this may be deeper than most damage sites on the Shuttle in flight, such a 
cavity would produce a stronger heating signal that is required to appropriately address the effect of the pressure 
gradient, and to develop an accurate assessment of the associated uncertainties. For both cavities, width is selected 
as W/H=3.6 to enhance the optical access needed for the global phosphor thermography used in this study. 

VI. Experiment Design 

This section outlines the experiment design process. First, the design philosophy is presented, along with 
physical constraints imposed on the process. This is followed by the development of the model test- surface 
geometries using two-dimensional CFD simulations, then a full three-dimensional CFD assessment that is used as a 
consistency check is presented. In the final section, the cavity design is generated, yielding the corresponding, 
available test space for these models. 

A. Design Philosophy and Constraints 

Two different nearly- constant pressure gradient test surfaces were developed to isolate the influence of pressure 
gradient on cavity heating augmentation. The first used a flat- surface geometry where dC p /ds~0.0, allowing 
consistency with existing data; the other provided a gradient near the target value of dC p /ds~-0. 15/ft (-0.013/inch). 
As these surfaces were being developed, the variation of the surface flow properties were mapped as a function of 
model test conditions and attitude (i.e. Pti, Tti, a, M e , and Re e ). Next, the cavity test location on the plate, X cav , was 
established downstream of the rapidly varying pressure gradients that occur as flow expands over the leading edge. 
This location was selected so that the desired scaled cavity with prescribed (H/8, W/H, L/H) yielded reasonable 
dimensions in geometric space (H, W, L) for measurement resolution. Because of the uniqueness of each model 
surface, there is only a single test point 
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(M e , Re e , H/8, W/H, L/H) Flat 

Surface Cavity (M e , Re e , H/8, W/H, L/H) GradieiU 

Surface Cavity 

where the cavity entry conditions were matched. Also, it is important to note that the available test domain for each 
model was developed naturally as a part of the design process. 

Other design constraints were considered in addition to matching cavity entry conditions. First, a blunted nose 
model was required to reduce the edge Mach numbers to target values at smaller angles of attack, thus minimizing 
the possibility of partial or complete tunnel-flow blockage (Dye, Everhart, and Hodge 1992). Effect of nose radius 
(Rn) was assessed as a part of the baseline flat plate design process. Second, the model width was made large 
enough to remove possible interference of side-edge vortices on the model test region and to minimize spanwise 
flow-gradient effects across the model. According to Micol (1995), the 31 -Inch Mach 10 Tunnel test core is safely 
10 inches at the smallest freestream Reynolds number. Therefore, a model width of 10 inches was selected. Finally, 
with proper mounting, models as long as 30 inches can be injected through the tunnel sidewall. However, there are 
low-pressure limitations on flow expansion for negative pressure gradient models that further limit the length of the 
model (the pressure can not be reduced below zero). Model length was defined during the design as nominally 20 
inches. 

B. Two-Dimensional Surface Modeling 

Two-dimensional (2D) viscous modeling of the test surfaces was performed using the LAURA computational 
fluid dynamics code (Gnoffo, et al., 1989, 1990). Solutions for the flat surface were computed for leading edge nose 
radius values (Rn) of 0.0625, 0.1250, and 0.1875 inches for a range of stream conditions covering tunnel Reoo = 0.5 
xlO 6 , 1.0 xlO 6 , 1.8 xlO 6 , and 2.0x1 0 6 at anticipated model angles of attack, a, ranging from -25 to 10 degrees. Most 
of the expected initial large change in the pressure was completed by the 5 -inch station for all conditions, 
particularly for the expanding flow conditions that occur around the test surface for positive angles of attack, though 
some slight gradient effects remain as the boundary layer continued to grow on the plate. The smallest nose radius 
was not selected because of fabrication concerns. Overall, relatively small differences were observed when 
comparing the effects of the 0.1250-inch nose radius to the 0.1875-inch nose radius on the boundary layer edge 
conditions. For consistency with other experiments, the 0.1250-inch nose was selected. Sample surface flow 
computations for the flat surface are present in Fig. 3. 
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Figure 3.- Flat plate surface flow properties for 
Re m =lxl0 6 . 



The challenge of designing an expansion surface was to maintain a linear longitudinal pressure distribution along 
the extent of the model, yielding a constant gradient, dC p /dx of the required magnitude over the cavity test region. 
Initially, power law surfaces of the form (z/z re f)=(x/L ref ) n were assumed with constant n values. However, the 
pressure gradient over these surfaces was not constant, though the surface did yield a nearly constant Mach number 
gradient. Subsequently, a reasonably constant pressure gradient was obtained by allowing n to vary parabolically 
along the surface according to n=2.0+1.5(x/L re f) 2 . Values of 3 inches and 26.55 inches were used for z re f and L re f, 
respectively. Sample computational results for the final geometry are presented in Figure 4. Solid colored lines are 
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used to represent the numerical solutions for different angles of attack with the exception of those on the pressure 
versus x subplot, where a consistent line-color-matched set of symbols is used. The solid lines through the symbols 


on the pressure versus x subplot represent a linear least 
squares fit from x=6 to x=19 inches. As can be observed, 
very satisfactory results for the linearity are achieved. 

The pressure gradient plotted versus angle of attack is 
present in Figure 5 for both the flat and expansion 
surfaces. Notably, for the expansion model the gradient 
is insensitive to the stream Reynolds number; but is 
highly dependent on the angle of attack. For the flat 
model, the gradient variation is very small with respect to 
a and is insensitive to Re (as with the expansion model). 
Prescribing a specific value of dC p /dx at a given set of 
cavity entry conditions (i.e. M e and Re 0 ) ties the results 
to a specific location on the vehicle surface, which goes 
beyond the binary yes/no question being addressed here. 
To achieve such a set of conditions would require 
additional iterations of the test surface geometry. 



Figure 5.- Pressure gradient variation with a for 
flat and expansion models. 


C. Three-Dimensional Flow Assessment 

The center of the plates are expected to be 
predominantly 2D flow in the streamwise direction, but 
three-dimensional (3D) flow relief is known to exist 
near the plate edges and the spanwise extent of this 
non-uniformity must be established. Ideally, 
instrumented pressure models would be tested to 
directly measure the baseline pressure gradients to 
check the validity of the two-dimensional design 
assumptions; however, resource limitations prevented 
this approach. An alternative approach was to perform 
a fully three-dimensional numerical assessment using 
the LAURA CFD code. 

To perform the 3D analysis, the support strut, which 
connects the leeward surface of the test article to the 
wind tunnel, was not modeled, but the expansion- 
surface plate and the mounting plate between the test 
article and the strut were retained. This produces a 
curved surface similar to a super- critical airfoil shape, 
shown in Figure 6. Additionally illustrated is the 0-0 
grid topology used to map a 4.5 million point, 120 
block decomposition on the plate to accommodate 
gradients around the corners of the plate and provide adequate resolution of the flow in high gradient regions. The 
structured grid was dimensioned to provide the necessary doubling in grid dimensionality to enable grid 
convergence studies required for the assessment of solution accuracy. Three grid resolutions were used, termed 
coarse, medium, and fine, where each was eight times as large as the previously smaller mesh. Grid convergence for 
this study was obtained when the heat transfer, or heating, at the windside surface did not change by more than 5% 
in the region of interest. Solution convergence was obtained when the L 2 norm of the residuals was at least 6 orders 
of magnitude less than where the solution started and the change in heating was less than 1% between 2000 
iterations with sufficiently large relaxation rates for both viscous and inviscid terms. Utilizing these conditions, the 
grid-converged state of the centerline pressure distributions were compared to the 2D flow analyses in Figure 7. 



Figure 6.- Surface grid topology used to model the 
constant pressure-gradient expansion plate. 
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Figure 8.- Surface heat transfer coefficient for 
a=16° and Reoo=1.3xl0 6 /ft. 
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For design of the plate, the 2D and 3D flow 
simulations are in good agreement, each yielding a 
nearly constant pressure gradient along the streamwise 
direction of the flow. Referring to the surface heat 
transfer coefficient distributions in Figure 8, the cavities 
will be located in the white box region on the plate 
centerline. (Note that the longitudinal axis is reversed in 
this figure because of a coordinate system difference in the CFD formulation, yielding the leading edge at x=20 
inches and the trailing edge at approximately x=l inch.) According to the constant heating contours, 2D effects are 
maintained a minimum of half an inch from the centerline with some expected spanwise 3D degradation beyond. 
Complete geometry computations with and without a centerline symmetry condition were evaluated and indicated 
no significant change from the half plate solutions used for the present example. Therefore, use of the plate for 2D 
flow experimental work should pose no issue with respect to the flow field and streamwise pressure gradient. 


D. Test Condition Selection and Cavity Design 

Next, the location of the leading edge of the cavities needed to be determined along with the final cavity 
dimensions. The 8-inch station was selected as the cavity leading edge position based on the following factors: 1) 
the constant pressure gradient location begins at approximately 5 inches (see Figs. 3 and 4), 2) thicker boundary 
layers allow larger cavities, 3) the target- condition requirement for the magnitude and range of Re e , and 4) the 
downstream narrowing of the plate side-edge vortices. Test conditions at this location were interpolated from the 
CFD surface flow distributions, and they are presented in terms of performance charts for the flat and expansion 
surfaces in Fig. 9 and Fig. 10, respectively. In these figures, M e , Re 0 , 8, and 0 are plotted versus P T i for various 
angles of attack. The effect of T T i for the range of conditions considered was second order; therefore, it was 
neglected when determining the boundary layer edge properties. For each test surface, three nominal edge Mach 
number conditions (2.25, 2.50, and 2.85) covering the target range were specified. At each edge Mach number 
condition, the available Re e space was determined from the curves, along with the a required to achieve these entry 
conditions; then, the corresponding 8 and 0 were determined to help scale the cavities. Because flow develops 
differently over an expansion plate from that developing over flat plate, the tunnel conditions and a were unique to 
each model. Likewise the actual geometric size of each cavity was unique, because the value of 8 at x=8 was 
different for each model at the prescribed entry conditions. The boundary layer edge conditions were re-evaluated 
using the actual test conditions, and these conditions were used for data reduction and post-test analysis. 

Surface flow conditions derived from the performance charts suggested that target edge conditions of M e =2.5 
and Re 0 =22O could be readily achieved on both the model surfaces and that a centered range of off-design 
conditions would be available. The resulting boundary layer thickness at this condition was used as the scaling 
factor for the cavity depth, H/8. The as-built cavity dimensions are presented in Table 2. The original intent was to 
develop baseline models to acquire heating, pressure, and oil flow information to characterize the undisturbed 
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environment and a measured pressure gradient; however, only the baseline heating model was built and tested due to 
resource limitations. 
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Figure 9.- Flat plate performance chart. 


Figure 10.- Expansion model performance 
chart. 


Table 2.- As-built cavity geometries. 


Model 

H 

(inch) 

(inch) L(inch) 

W/H 

L/H 


CLOSED - Long Cavity 



Expansion 

0.1443 

0.5268 2.9011 

3.65 

20.1 

Flat 

0.2095 

0.7474 4.1622 

3.57 

19.9 


OPEN - Short Cavity 



Expansion 

0.1462 

0.5218 1.0445 

3.57 

7.1 

Flat 

0.2095 

0.7496 1.4969 

3.58 

7.1 


Experimental Method 

The experimental methods are presented here, beginning with a description of the test facility and nominal tunnel 
flow conditions. This is followed by a description of the model fabrication process, a discussion of the phosphor 
coating used for making the heating measurements, and the model mounting. An overview of the phosphor data 
system used to acquire the global heating measurements is next presented, followed by the post-test data mapping 
and scaling. The test data are presented in tabular form, first as a summary of the cavity and run parameters for each 
model, and then as global images of the surface heating. Finally, an uncertainty analysis is developed and sample 
results are presented in the form of global images of the bias and total uncertainty. 

A. Facility 

The 31 -Inch Mach 10 Air Tunnel was selected as the test facility because the higher surface temperature rise 
during a run would provide lower uncertainties in heat transfer measurements (Merski 2001). This is a blowdown 
facility in which heated, filtered air is used as the test gas. The tunnel has a square, contoured nozzle, which opens 
into a 31 -inch square test section. Models are supported on a hydraulically operated, sidewall-mounted injection 
system that can transfer a model from the sealed model box to the tunnel center-line in less than 0.6 seconds. Tunnel 
run times of approximately 60 seconds can be achieved, but typical heating studies require only a few seconds. The 
nominal reservoir conditions are stagnation pressures of 150 psi to 1450 psi at stagnation temperatures of about 
1825°R. These reservoir conditions very nearly produce perfect gas (y = 1.4) freestream flows with a Mach number 

of approximately 10 and Reynolds numbers of 0.2x1 0^/ft to 2.2x1 0^/ft. Nominal flow conditions are provided in 
Table 3 and a detailed description of this facility is presented by Micol (1995). 
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Table 3.- Nominal flow conditions for the 31-Inch Mach 10 Air Tunnel. 


P T 1 

psi 

T T i 

°R 

P„xl0 2 

psi 

T^ 

°R 

q~ 

psi 

ft/s 

M. 

R*,xl0‘ 6 

ft' 1 

Pt,2 

psi 

Test Core 
in. 

348 

1800 

0.992 

95.2 

0.650 

4628 

9.67 

0.568 

1.205 

10x10 

723 

1825 

1.867 

94.2 

1.259 

4670 

9.81 

1.104 

2.334 

12x12 

1452 

1800 

3.509 

90.7 

2.430 

4643 

9.94 

2.240 

4.503 

12x12 


B. Models 

The flat plate and expansion plate models used in this experiment were nominally 10-inches wide by 20-inches 
long. All models have a 0.125-inch radius nose with a 15-degree wedge on the lower surface. The general 
manufacturing process for ceramic test articles is 
described in Buck, et al. (1993). Specific changes to 
this process for a previous, unpublished cavity heating 
experiment are described in Buck, et al. (2006). In this 
later manufacturing method, a rapid prototyping 
stereolithography (SLA) system at NASA Langley 
Research Center was used to make a resin mold 
pattern for casting cavities in flat 4-inch wide by 18- 
inch long ceramic insert models. These inserts were 
then installed in a 10-inch wide by 28-inch long steel 
flat-plate model for testing. The present experiment 
further improves this manufacturing process by 
eliminating the insert requirement. Here, the entire 10- 
inch wide by 20-inch long ceramic upper surface is 
cast as a single surface, including the nose, and it is 
then bonded to an aluminum backing plate. This 
improvement removes the step and gap mounting 
problems experienced during installation for testing 
with the previous method. The backing plate is at 0° 
incidence relative to tunnel centerline for the flat plate, and it is offset at 8 U nose-up incidence for the expansion 
model. The cavities were machined into the plates, yielding sharp, well-defined edges and surfaces. These models 
were then coated with a mixture of phosphors suspended in a silica-based colloidal binder and sent to quality 
assurance for measurement and application of small circular locating markers, known as fiducial marks. These 
fiducial marks aid in data reduction and model orientation, and they can be seen in the run images as dark dots. Even 
though fiducial marks do not influence the flow over the model surface, they will influence the local heating 
measurement because the ink used in the process changes the calibration of the phosphor coating at their application 
point. Their influence on the data, notably in line cuts passing through them, is to produce a very short, very sharp- 
peaked sinusoidal signature covering the diameter of the marking. When mounted in the tunnel, the aluminum 
backing plate was attached to a wedge block surface (0° for the flat plate model and 8° for the expansion model) 
used to align the model to zero incidence. Finally, this entire assembly was attached from the leeward side to an 
existing model support strut and mounted to the tunnel sidewall injection plate. An offset angle wedge is placed 
under the strut base when required to extend the angle testing range of the model. A picture of the model mounted in 
the 31 -Inch Mach 10 Tunnel with the phosphor thermography camera and lighting system is given in Fig. 11. 



Figure 11.- Test model installed in 31-Inch Mach 10 
Air Tunnel with TGP lighting and camera system. 


C. Thermographic Phosphor Measurement System 

Global surface heating distributions were calculated using the two-color, relative-intensity, phosphor- 
thermography aeroheating measurement method (Merski 1999). This is the standard method for obtaining 
aeroheating data in NASA Langley’s hypersonic wind tunnels, and it can be used to identify the surface heating 
effects of complex three-dimensional flow phenomena, which are difficult to examine using conventional discrete- 
sensor methods. With this method, ceramic wind tunnel models are coated with phosphor crystals that fluoresce in 
the red and green regions of the visible light spectrum when illuminated by ultraviolet (UV) light. During a wind 
tunnel run, the phosphor-coated model is exposed to the heated flow of the tunnel, and the resulting changes in 
fluorescence intensity of the model are recorded and digitized through a 640 by 480 resolution color CCD (charge 
coupled device) camera and a state-of-the-art video acquisition system. The fluorescence intensity is dependent on 
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both the intensity of the incident UV light and the local model surface temperature. The UV dependence is removed 
by taking the ratio of the green to red intensity images, from which surface temperature distributions are determined 
through prior calibrations. Images are acquired before the wind tunnel run and after injection of the model to the 
tunnel centerline during a run. 

D. Data Reduction 

Global mappings of the surface temperature obtained with the thermographic phosphor (TGP) data acquisition 
system were reduced to surface heating distributions using the IHEAT data reduction software (Merski 1999). In this 
method, phosphor images are acquired shortly after injection of the model to the tunnel centerline, which typically 
requires less than one second. With this software, temperature mappings acquired through the run are used to obtain 
heat transfer mappings. This is done by applying one-dimensional, semi-infinite-solid heat conduction theory 
assuming a constant heat- transfer coefficient, and by making empirical corrections to account for temperature 
changes in model substrate thermal properties. The results are presented in terms of a non-dimensional heat transfer 
coefficient ratio, h/h re f, where h re f is generally taken as h FR , the theoretical heating computed with the Fay-Riddell 
(1958) theory. For this experiment, h FR was computed using a 0.125-inch radius sphere (the test model nose radius) 
and a reference stagnation temperature of 540°R. At this point, IHEAT was used to extract preliminary analysis data 
along lines from the two-dimensional heating images. Also, these images were mapped onto a three-dimensional 
representation of the test article using IHEAT ’s MAP3D photogrammetry program for further global analysis. 

Data for this experiment were translated such that the leading edge of the cavity is at x=0. The (x, y, z) data were 
then scaled to (x/H, y/H, z/8) using the measured cavity depth, H, and the computed boundary layer thickness, 8. 
This yields a scaled-geometry cavity of depth (H/8), width (W/H), and length (L/H). Note again that the boundary 
layer edge conditions are determined from the LAURA computations for matching tunnel and model conditions for 
each run. For the contoured expansion models where z is not zero, a baseline longitudinal profile outside of the 
cavity is subtracted from the contour containing a cavity, yielding an incremental data surface where z=0 
everywhere except in the cavity. 

The impact of the cavity on the local undisturbed environment was assessed by converting the local heating data 
to heating augmentation or bump factor (BF) format by normalizing (hi oca i/h ref ) by a reference- location average 
heating (h a vg/h re f), yielding BF=(hi oca i/h avg ). In this format, BF=1 is the nominal undisturbed condition for a flat plate; 
however, in an expanding flow BF will decrease slightly down the plate as the heating drops. For analysis 
consistency, (h avg /h re f) is computed upstream of the cavity over an area two cavity depths long by one cavity width 
wide, and this region is located at (-3 < x/H < -1) by (-0.5 W/H < y/H < 0.5 W/H). 

E. Presentation of Test Data 

For completeness, the as-run test conditions for the as-built cavity geometries are presented in Table 4 through 
Table 7. Each table represents the block of runs dedicated to a single configuration, providing the computed cavity 
entry conditions M e , Re 0 , Re 0 /M e , 8, and 0, as well as the cavity depth, width, and length normalized by the 
boundary layer thickness. The pressure gradient, dC p /dx, is also provided. Each configuration contains a randomly 
selected repeat run as a consistency check on the data acquisition and reduction processes. The corresponding 
mapped and scaled bump factor images are presented in Table 8 through Table 11 with the corresponding run 
number identified on each image. Following the usual convention, flow is from left to right. Note on some runs that 
the camera field of view truncates the image, introducing a cold arc- shaped region at the upstream top corner edge 
and a cold strip on the downstream edge spanning the model width. Examples of this truncation are visible in runs 
23 through 25 and runs 55 through 66. 

The images shown in Tables 8-11 provide a wealth of both local and global information, representing a 
significant analysis task unto itself, so only a few observational remarks are provided here. Flow entering the long, 
closed cavities (runs 23 through 44) is particularly noticeable in the augmented heating levels on the cavity floor. 
These cavities also exhibit enhanced vortical streaking and wiping that increases the local surface shearing, leading 
to large heating values on the cavity end wall and further downstream. In contrast and as expected, the short, open 
cavities (runs 45 through 66) are much cooler on the cavity floor. Large heating augmentation is only dominant on 
the downstream end wall lip. Though vortical streaking is still visible aft of the cavity, it decays rapidly. All of these 
observations are in accordance with published literature. Particularly noteworthy observations are that the location 
of peak heating may occur in the corners; depending on the flow conditions, peak heating will move laterally; and, 
discrete sensors placed on the centerline would in all likelihood have missed the peak heating values that are readily 
observed when using TGP methods. 
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Table 4 - Expansion Plate Model - L/H=20, W/H=3.6 Cavity. 


Run 

Model 

M e 

Re e 

Re e /M e 

8 (inch) 

0 (inch) 

dCp/dx 

(1/ft) 

H/8 

W/8 

L/8 

23 

Exp 

2.848 

151 

53 

0.2340 

0.0275 

-0.071 

0.62 

2.3 

12.4 

24 

Exp 

2.872 

171 

60 

0.2070 

0.0242 

-0.063 

0.70 

2.5 

14.0 

26 

Exp 

2.897 

211 

73 

0.1653 

0.0193 

-0.055 

0.87 

3.2 

17.6 

27 

Exp 

2.910 

190 

65 

0.1844 

0.0217 

-0.057 

0.78 

2.9 

15.7 

28 

Exp 

2.377 

196 

82 

0.1606 

0.0211 

-0.159 

0.90 

3.3 

18.1 

29 

Exp 

2.441 

213 

87 

0.1402 

0.0186 

-0.131 

1.03 

3.8 

20.7 

30 

Exp 

2.487 

234 

94 

0.1283 

0.0166 

-0.110 

1.12 

4.1 

22.6 

31 

Exp 

2.499 

254 

102 

0.1173 

0.0152 

-0.110 

1.23 

4.5 

24.7 

32 

Exp 

2.309 

214 

93 

0.1409 

0.0194 

-0.183 

1.02 

3.7 

20.6 

33 

Exp 

2.353 

234 

100 

0.1245 

0.0167 

-0.144 

1.16 

4.2 

23.3 

34 

Exp 

2.414 

254 

105 

0.1154 

0.0152 

-0.124 

1.25 

4.6 

25.1 



Table 5. 

- Flat Plate Model - 

L/H=20, W/H=3.6 Cavity 

. 



Run 

Model 

M e 

Re e 

Re e /M e 

8 (inch) 

0 (inch) 

dCp/dx 

(1/ft) 

H/8 

W/8 

L/8 

35 

Flat 

2.825 

186 

66 

0.2166 

0.0259 

-0.004 

0.97 

3.4 

19.2 

36 

Flat 

2.847 

204 

72 

0.1961 

0.0234 

-0.004 

1.07 

3.8 

21.2 

37 

Flat 

2.873 

224 

78 

0.1811 

0.0214 

-0.005 

1.16 

4.1 

23.0 

38 

Flat 

2.900 

244 

84 

0.1691 

0.0197 

-0.005 

1.24 

4.4 

24.6 

39 

Flat 

2.407 

220 

91 

0.1634 

0.0219 

0.000 

1.28 

4.6 

25.5 

40 

Flat 

2.461 

241 

98 

0.1477 

0.0194 

0.000 

1.42 

5.1 

28.2 

41 

Flat 

2.478 

260 

105 

0.1347 

0.0175 

0.000 

1.55 

5.5 

30.9 

42 

Flat 

2.170 

278 

128 

0.1242 

0.0178 

-0.001 

1.69 

6.0 

33.5 

43 

Flat 

2.205 

301 

136 

0.1107 

0.0155 

0.000 

1.89 

6.8 

37.6 

44 

Flat 

2.461 

241 

98 

0.1477 

0.0194 

0.000 

1.42 

5.1 

28.2 


Table 6.- Expansion Plate Model - L/H=7.2, W/H=3.6 Cavity. 


Run 

Model 

M e 

Re e 

Re e /M e 

8 (inch) 

0 (inch) 

dCp/dx 

(1/ft) 

H/8 

W/8 

L/8 

55 

Exp 

2.848 

151 

53 

0.2339 

0.0437 

-0.071 

0.62 

2.2 

4.5 

56 

Exp 

2.580 

191 

74 

0.1648 

0.0210 

-0.063 

0.89 

3.2 

6.3 

57 

Exp 

2.910 

190 

65 

0.1844 

0.0217 

-0.057 

0.79 

2.8 

5.7 

58 

Exp 

2.897 

211 

73 

0.1653 

0.0193 

-0.055 

0.88 

3.2 

6.3 

59 

Exp 

2.377 

196 

82 

0.1606 

0.0211 

-0.159 

0.91 

3.2 

6.5 

62 

Exp 

2.442 

213 

87 

0.1406 

0.0187 

-0.131 

1.04 

3.7 

7.4 

63 

Exp 

2.487 

233 

94 

0.1287 

0.0167 

-0.110 

1.14 

4.1 

8.1 

61 

Exp 

2.499 

254 

102 

0.1173 

0.0152 

-0.110 

1.25 

4.4 

8.9 

64 

Exp 

2.309 

214 

93 

0.1409 

0.0194 

-0.183 

1.04 

3.7 

7.4 

65 

Exp 

2.353 

234 

100 

0.1245 

0.0167 

-0.144 

1.17 

4.2 

8.4 

66 

Exp 

2.414 

254 

105 

0.1154 

0.0152 

-0.124 

1.27 

4.5 

9.1 

67 

Exp 

2.442 

213 

87 

0.1406 

0.0187 

-0.131 

1.04 

3.7 

7.4 



Table 7. 

- Flat Plate Model - L/H=7.2, W/H=3.6 Cavity 

. 



Run 

Model 

M e 

Re e 

Re e /M e 

8 (inch) 

0 (inch) 

dCp/dx 

(1/ft) 

H/8 

W/8 

L/8 

45 

Flat 

2.825 

186 

66 

0.2167 

0.0278 

-0.004 

0.97 

3.5 

6.9 

46 

Flat 

2.847 

204 

72 

0.1961 

0.0234 

-0.004 

1.07 

3.8 

7.6 

47 

Flat 

2.873 

224 

78 

0.1810 

0.0213 

-0.005 

1.16 

4.1 

8.3 

48 

Flat 

2.900 

244 

84 

0.1691 

0.0197 

-0.005 

1.24 

4.4 

8.9 

49 

Flat 

2.407 

220 

91 

0.1633 

0.0219 

0.000 

1.28 

4.6 

9.2 

50 

Flat 

2.461 

241 

98 

0.1477 

0.0194 

0.000 

1.42 

5.1 

10.1 

51 

Flat 

2.477 

260 

105 

0.1347 

0.0175 

0.000 

1.56 

5.6 

11.1 

52 

Flat 

2.170 

278 

128 

0.1241 

0.0178 

-0.001 

1.69 

6.0 

12.1 

53 

Flat 

2.205 

301 

136 

0.1107 

0.0155 

0.000 

1.89 

6.8 

13.5 

54 

Flat 

2.825 

186 

66 

0.2166 

0.0259 

-0.004 

0.97 

3.5 

6.9 


14 of 28 

AIAA Paper 2006-0185 


































Table 8.- Expansion Plate Bump Factor Images - L/H=20, W/H=3.6 
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Table 8.- Concluded. 



Test 406 Run 30 


i 


Bump 

10.00 

7.50 

5.00 

4.00 

3.50 

3.00 

2.50 

2.00 

1.75 

1.50 
1.25 
1.00 
0.75 
0.50 
0.25 
0.00 



Test 406 Run 31 


I 


Bump 

10.00 

7.50 

5.00 

4.00 

3.50 

3.00 

2.50 

2.00 

1.75 

1.50 
1.25 
1.00 
0.75 
0.50 
0.25 
0.00 



Test 406 Run 32 


i 


Bump 

10.00 

7.50 

5.00 

4.00 

3.50 

3.00 

2.50 

2.00 
1.75 

1.50 
1.25 
1.00 
0.75 
0.50 
0.25 
0.00 



Test 406 Run 33 


i 


Bump 

10.00 

7.50 

5.00 

4.00 

3.50 

3.00 

2.50 

2.00 

1.75 

1.50 
1.25 
1.00 
0.75 
0.50 
0.25 
0.00 



Test 406 Run 34 


i 


Bump 

10.00 

7.50 

5.00 

4.00 

3.50 

3.00 

2.50 

2.00 

1.75 

1.50 
1.25 
1.00 
0.75 
0.50 
0.25 
0.00 
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Table 9.- Flat Plate Bump Factor Images - L/H=20, W/H=3.6. 







I 


Bump 

10.00 

7.50 

5.00 

4.00 

3.50 

3.00 

2.50 

2.00 

1.75 

1.50 
1.25 
1.00 
0.75 
0.50 
0.25 
0.00 





Bump 

10.00 

7.50 

5.00 

4.00 

3.50 

3.00 

2.50 

2.00 

1.75 

1.50 
1.25 
1.00 
0.75 
0.50 
0.25 
0.00 
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Table 9.- Concluded. 



Table 10.- Expansion Plate Bump Factor Images - L/H=7.2, W/H=3.6. 
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Table 10.- Continued, 
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Table 10.- Concluded. 



Table 11.- Flat Plate Bump Factor Images - L/H=7.2, W/H=3.6. 
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Table 11.- Concluded, 








l 


Bump 

10.00 

7.50 

5.00 

4.00 

3.50 

3.00 

2.50 

2.00 

1.75 

1.50 
1.25 
1.00 
0.75 
0.50 
0.25 
0.00 


I 


Bump 

10.00 

7.50 

5.00 

4.00 

3.50 

3.00 

2.50 

2.00 

1.75 

1.50 
1.25 
1.00 
0.75 
0.50 
0.25 
0.00 


I 


Bump 

10.00 

7.50 

5.00 

4.00 

3.50 

3.00 

2.50 

2.00 

1.75 

1.50 
1.25 
1.00 
0.75 
0.50 
0.25 
0.00 


Bump 



I .JU 

1.25 


1.00 

0.75 

0.50 

0.25 

0.00 


21 of 28 

AIAA Paper 2006-0185 


F. Uncertainty Analysis 

Merski published the uncertainties for the phosphor thermography method as a whole in 1999. A more 
quantitative understanding of the uncertainty levels for this experiment were obtained by adding phosphor heat- 
transfer bump-factor uncertainty algorithms to the MAP3D code, resulting in global uncertainty surface maps 
similar to the three-dimensional heat transfer maps. The process is outlined below. 

The solution of the heat conduction equation (i.e. the data reduction equation) used in the IHEAT code is given 
by 


0 ,„ 


(T w /iJ-T ini 


= 1 - e'erfcA 


where 


Q(x, y,t) = T (x, y, t ) - T inil (x, y,0) and 


A KhJTjft 
P 


Also, h is the convective heat transfer coefficient, h aw is the adiabatic wall enthalpy, h w is the wall enthalpy, T w is 
the corresponding wall temperature, and t is the effective time of data acquisition. /? is the thermal product of the 
substrate materials and it is the square root of the product of the material density, specific heat, and thermal 
conductivity. Bias and precision uncertainty values developed at the 95 percent confidence level for each parameter 
appearing in the data reduction equation are given in Table 12. 


Table 12.- Variable uncertainty values. 


Uncertainty Type 

Bias 

Uncertainty Value 

Precision 

Uncertainty Value 

Initial model wall temperature 

1.43 °C 

1.0 °C 

Model run wall temperature at reference location 

1.93 °C 

n/a 

Model run wall temperature 

3.36 °C 

1.0 °C 

Effective time 

0.02 s 

0.05 s 

Thermal properties 

5.9% 

n/a 


Each of these individual parameter uncertainties was individually inserted into the data reduction equation to 
determine a component uncertainty in heat transfer coefficient. The bias uncertainty, B , and the precision 
uncertainty, P , for the heat transfer coefficient were determined by obtaining the root- sum-square (RSS) from each 
of the component uncertainties using 



( 


1/2 

( 


B = 

5 >/ 

and P = 

2 >/ 


v J 

) 


^ J 

J 


where j is the number of component uncertainties. Total uncertainties are obtained by taking the RSS of the bias 

and precision uncertainties. Bias and total uncertainties of the heat transfer coefficients are determined at every pixel 
point imaged on the model. 

Bump factors extracted from the mappings were used in two ways: 1) as line cuts when selecting specific data, 
and 2) as regions of interest over which all of the data were averaged. In the first case, the total uncertainties are 
applicable and in the latter, bias uncertainties are applicable. Precision uncertainties are removed during averaging, 
since they are primarily due to random pixel scatter. Accordingly, the average heat transfer coefficient h avg used to 
compute bump factors was assumed to have only a bias uncertainty. Typically, total uncertainties are determined 
using the quotient rule by adding the total uncertainty in hi oca i to the bias uncertainty in h avg . For nominally flat 
models of this type, very low heating rates are obtained during a run and the resulting uncertainties were initially 
very high. Yet, comparisons of line cuts were more consistent than the uncertainty analysis seemed to suggest. 

While conservatism is important, excessive conservatism is undesirable. Therefore, the phosphor temperature 
lookup table data were re-examined and it was determined that since there was minimal variation of incident UV 
intensities on the models during the tunnel runs, and because of the low heating measurements, the data were 
typically confined to one very small segment of the temperature range. Thus, the bias uncertainties for the data used 
in calculating hi oca i and h avg were uni-directional and determinate. Therefore, the bias uncertainties in hi oca i had to be 
subtracted from the bias uncertainties in h avg via the quotient rule, instead of added, when computing bias 
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uncertainties for the augmentation factors. Similarly, total uncertainties in augmentation factors were then obtained 
simply by taking the RSS of the bias and precision uncertainties in h avg and by subtracting the bias uncertainties in 

hlocal* 

Representative bias uncertainties for Run 63, an open cavity (L/H=7.2, H/8=1.14) on the expansion plate, are 
shown in Figure 12a, corresponding total uncertainties are shown in Figure 12b. This cavity represents one of the 
more challenging measurement configurations from an accuracy standpoint because the flow does not directly enter 
the cavity, resulting in a very small floor temperature rise relative to the pre-run condition. Here, total and bias 
uncertainties are typically less than 10% over the plate surface. In the cavity where the lowest temperatures are 
experienced, the bias uncertainties are about 10%, while the total uncertainties approach 20%. 



a) Bias uncertainties (% Uncertainty/100). 


b) Total uncertainties (% Uncertainty/100). 


Figure 12.- Expansion model bump factor uncertainties for Run 63 - L/H=7.2, H/S=1.14. 


VII. Results 

The objective of the present paper is to address whether pressure gradient has an impact on the augmentation of 
cavity-induced heating. The results presented here represent a preliminary evaluation of the data based solely a 
qualitative examination of the centerline BF distributions that are presented in Figure 13 and Figure 14. In each case, 
the scaled and normalized distribution data are presented with the cavity entrance at (x-x cav )/H=0. The cavity end 
wall occurs nominally at (x-x cav )/H=7.2 for the short (open) cavities and (x-x cav )/H=20 for the long (closed) cavities. 
The sharp spikes in the data upstream and downstream of the cavity are the disturbances caused by the previously 
mentioned fiducial markings on the model. Normally, these fiducial marks would be filtered out of the data, but they 
have been left in to indicate the quality of the data alignment that occurs during translation and scaling processes. 

Considering the short cavity cases first (Fig. 13), flow remains laminar for all the runs with the expansion model; 
however, runs 52 and 53 with the flat plate model appear non-laminar upon exiting the cavity. The favorable 
pressure gradient provided by the expansion plate offers some relief to the susceptibility of the short cavity to 
boundary layer transition. When transitional flow is observed downstream a different, elevated heating profile 
occurs inside the cavity, most likely because of a more energetic shear layer over the cavity, increasing the cavity 
mass interchange with the hotter main stream. Not unexpectedly, little difference is seen in the distributions when 
these cases are excluded, because the flow physics implies that the external stream skips the cavity. 

After first excluding transitional/turbulent runs 41, 42, and 43 occurring on the flat plate, the long cavity 
distributions presented in Figure 14 yield a different conclusion. For long (closed) cavities, the flow will enter the 
cavity. Since the pressure in the expansion model is dropping below that of the flat plate model, it is anticipated that 
the turning of the flow into the cavity will be greater and that vortical interactions will be increased. Figure 14a 
supports this conclusion based on the augmented heating levels on the floor of the cavity compared to the zero- 
pressure gradient cavity shown in Figure 14b. Similarly downstream of the expansion cavity, the increased vorticity 
in the flow elevates the shearing at the surface and augments the centerline heating above that experienced by the 
zero-gradient cavity. However, a fair comparison requires further global image analysis to establish whether the 
flow exiting the long expansion cavities has remained laminar. 
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13a- Expansion model. 13b) Flat model. 

Figure 13.- Centerline bump factor distributions for L/H=7.2 models. 




14a) Expansion model. 14b) Flat plate model. 

Figure 14.- Centerline bump factor distributions for L/H=20 models. 


A more quantitative assessment can be obtained by 
averaging the floor data between 0.5 < (x-x cav )/H < 6 for 
the short cavities and between 0.5 < (x-x cav )/H < 18 for 
the long cavities to create BF avg . These limits are 
selected to remove the upstream and downstream end 
wall effects on the data. The variation of BF avg with 
pressure gradient is presented in Figure 15. The circles 
represent the long cavity; squares represent the short 
cavity; and, filled symbols represent the nominally zero 
pressure gradient conditions on the flat plate model. 
Considering first the short cavities (squares), BF avg , the 
floor heating augmentation, is roughly 20% of the 
nominal undisturbed surface heating, irregardless of the 
pressure gradient. The long cavities (circles), on the 
other hand, have a distinct pressure gradient effect with 
the average heating on the cavity floor being about 30% 
of the reference for low-pressure gradients and 



Figure 15.- Variation of average centerline floor 
bump factor with pressure gradient. 
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Figure 16.- Variation of maximum centerline bump 
factor with pressure gradient. 


increasing to nearly 50% as the pressure gradient increases. Note that this increase is above the previously cited 
uncertainty level. 

Figure 16 presents the variation of the maximum 
end wall bump factor, BF max , as an additional 
pressure-gradient assessment metric. All runs, 
including those presumed transitional/turbulent have 
been included. No clear discernable trend with 
pressure gradient is apparent here for either the short 
or long cavity, but a large jump in the magnitude from 
the 3-to-5 range to the 8-to-10 range does occur for the 
long-cavity zero-gradient test conditions where 
transitional/turbulent flow is presumed (runs 41, 42, 
and 43). Plotting the BF max variation with H/8 offers 
no additional clarity; however, if the variation is 
plotted with L/8 (Fig. 17), then the step change occurs 
for L/8>25. Though not complete as a correlation 
parameter, L/8 is a clear indicator of cavity- induced 
transition/turbulence onset for those cavities installed 
on the flat plate. An examination of the traditional 
boundary layer correlation parameter, Re 0 /M e , found 
in Table 4 and Table 5 also indicates that these 
represent the largest values obtained during the 
experiment, and according to unpublished, ongoing 
independent work, they are in the range where 
transition can be expected. Further analysis of the 
image data and correlation with additional data are 
required before any definitive expression can be made 
regarding the effect of pressure gradient on BF max . 

VIII. Conclusion 

An experiment to evaluate the possible influence 
of pressure gradient on the heating augmentation 
induced by rectangular cavities has been conducted in 
support of the Space Shuttle Return-to-F light 
Program. Two-dimensional viscous computational 
solutions via the LAURA code were used to design a 
rectangular-planform expansion-plate model, yielding 
a longitudinal gradient typical of what would be 
experienced on the chine region of the shuttle orbiter. 

Three-dimensional viscous solutions were used to evaluate the centerline pressure gradient as a check on the design 
process. Using the two-dimensional simulations of the expansion plate model and a corresponding flat plate model, 
predicted boundary layer edge conditions were used to design open (L/H=7.2, W/H=3.6, H/8=l.l) and closed 
(L/H=20, W/H=3.6, H/8=l.l) cavities for each model. Four ceramic cavity models (one for each cavity) were 
manufactured and the two-color phosphor thermography method was used to determine the global heating 
distribution on each. The test parameters and image data have been presented for completeness. Centerline 
distributions normalized to bump factor (or heating augmentation factor) format have been analyzed to assess the 
existence of a pressure gradient effect. For the short (open) cavity geometries, pressure gradient was found to have 
minimal effect, if any, on the average floor heating augmentation. In this case, all floor averages scattered around a 
nominal value of 0.2. For the long (closed) cavity geometries, pressure gradient increased the average heating 
augmentation factor by 50% on average above the zero gradient conditions, which was approximately 0.3. An 
assessment of the maximum augmentation on the end wall revealed no apparent effect for the short cavity; however, 
the analysis was inconclusive for the long cavity, requiring additional image and correlation analysis. 



L/8 


Figure 17. Variation of maximum centerline bump 
factor with L/8. 
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